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В данной статье рассматривается решение задачи диагностики дорожных покрытий на реконфигурируемой 
вычислительной системе с применением языка СОГАМО. Подробно рассмотрены особенности реализации 
вычислительно трудоемких частей задачи, представлены граф-схемы, приведены примеры текста программы 
на языке СОГАМО. 


С помощью электромагнитного зондирования удается оперативно получать струк- 
туру дорожного покрытия, определять толщины составляющих покрытие слоев. В зави- 
симости от используемой частоты антенны подход позволяет получать структуру 
дороги с различными разрешением и глубиной. Метод позволяет идентифицировать 
слои асфальта, бетона, песка, глины, пустоты, переувлажненные участки. 

Постановка задачи. В задаче диагностики дорожных покрытий можно выде- 
лить несколько последовательных этапов обработки, которые проходит каждая трасса, 
поступающая из георадара, заключающихся в решении параметрических математи- 
ческих задач, решении интегрального уравнения, применении методов обработки 
изображения и распознавания образов. Решением задачи является вертикальное се- 
чение дороги (под линией зондирования) в виде обработанного полутонового изобра- 
жения, а также выделенные в процессе обработки слои и их толщины. 

Анализ алгоритма задачи показал, что наибольшей вычислительной трудоем- 
костью обладают этапы вычитания прямого сигнала и обращения интегрального 
уравнения. Все остальные этапы и другие вычислительные блоки могут быть выпол- 
нены на обычном персональном компьютере. Таким образом, схема реализации за- 
дачи диагностики дорожных покрытий для реконфигурируемых вычислительных 
систем (РВС) состоит из двух вычислительно трудоемких модулей, каждый из кото- 
рых представляет собой отдельную программу на языке СОГАМО. 

Язык программирования СОГАМО представляет собой высокоуровневый язык, 
предназначенный для программирования многопроцессорной вычислительной систе- 
мы с реконфигурированной архитектурой. 

Язык СОГАМО позволяет описывать вычислительные алгоритмы для реализа- 
ции на реконфигурируемых вычислительных системах в пределах одной архитектуры 
различных конфигураций и дает возможность программисту максимально просто опи- 
сывать различные виды параллелизма в достаточно сжатом виде. 

Рассмотрим реализацию модуля, вычитающего прямой сигнал, на РВС. Каждая 
трасса предварительно нормируется с помощью умножения каждого отсчета трассы 
на нормировочный коэффициент. После этого осуществляется поиск оптимальных 
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значений сдвига и коэффициента вычитания методом перебора сдвига в заданном 
диапазоне с помощью вычисления кубического сплайна, интерполирующего задан- 
ные значения табличной функции в предварительно рассчитанных точках. В конце 
поиска на основе найденных оптимальных значений сдвига и коэффициента вычитания 
вычисляется прямой сигнал, который, в свою очередь, вычитается из нормированной 
трассы. 

На рис. 1 представлена граф-схема, осуществляющая перебор сдвига в диапазоне 
для нахождения оптимального значения сдвига и коэффициента для вычитания пря- 
мого сигнала. 
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Рисунок 1 — Граф-схема поиска оптимального сдвига и коэффициента 
для вычитания прямого сигнала методом перебора сдвига 


Граф-схема вычислений - это структура, отображающая информационную за- 
висимость вычислений для каждого элементного узла. Каждая вершина этой структуры 
является элементарным процессором (функциональным устройством), участвующим 
в решении задачи. 

На рис. 1 приняты следующие обозначения: 

— Г - генератор; 

— Р1: - исходные параметры для текущего шага перебора для текущей трассы; 
— $ — блок вычисления суммы элементов потока в заданном диапазоне; 

— БВСИ - блок вычисления сплайн-интерполяции; 

— 2 — блок возведения в квадрат; 

— 7 — блок задержки данных; 

— С; - константы; 

— БВЗ — блок выборки значений; 

— К — значение коэффициента для вычитания прямого сигнала для текущей трассы; 
— М - значение невязки для текущей трассы; 

— Р2; - выходные параметры для текущего шага перебора для текущей трассы. 

На каждом шаге перебора для каждой трассы вычисляется невязка М. Если это 
первый проход (изначально минимальная невязка отрицательна) или вычисленная 
невязка меньше минимальной, то (при условии, что еще не выполнены все шаги пе- 
ребора) в выходные параметры записываются текущее значение сдвига, вычислен- 
ный коэффициент вычитания К, вычисленная невязка М, иначе записываются данные 
с предыдущего шага. Коммутацию выходных значений обеспечивает блок выборки 
значений, который представляет собой набор мультиплексоров, управляемых сигна- 
лом 5е[. 

В задаче диагностики дорожных покрытий область оптимизации включает в себя 
60 отсчетов трассы. Дополнение области нулями до 64 элементов позволяет находить 
сумму этих элементов, используя схему, представленную на рис. 2. 
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Рисунок 2 — Граф-схема сложения 64 значений 


Блоки 2; задерживают данные на 1 тактов. На первом сумматоре происходит по- 
парное сложение 64 значений с получением 32 частичных сумм. Эти частичные суммы 
идут через такт, поэтому перед вторым сумматором необходимо увеличить задержку 
данных до двух. В результате получается 16 частичных сумм, но уже идущих каж- 
дый четвертый такт. Таким образом, увеличивая каждый раз задержку данных перед 
последующим сумматором вдвое, мы на выходе последнего сумматора получим сумму 
64 значений на 64-м такте. 

Блок вычисления сплайн-интерполяции позволяет находить значения кубичес- 
кого сплайна, используя табличную функцию в узлах сетки. Данный блок естественным 
образом реализуется на языке СОГАМО. Отметим, что сетка, табличная функция и 
коэффициенты кубического сплайна реализованы в виде блоков ПЗУ, так как в пре- 
делах обрабатываемых трасс остаются неизменными. 

Для организации плотного потока данных через вычислительную схему предла- 
гается осуществлять перебор не для каждой трассы в отдельности, а для всех трасс. 
То есть один кадр работы схемы будет представлять собой один шаг перебора для каж- 
дой трассы. Однако в данном случае количество кадров будет определяться относительно 
той трассы, которая имеет наибольшее количество повторений для перебора. 

Теперь рассмотрим реализацию модуля, решающего задачу обращения интеграль- 
ного уравнения на РВС. Данный этап обработки представляет собой операцию сверт- 
ки — над исходной трассой выполняется операция быстрого преобразования Фурье 
(БПФ), результат которого комплексно умножается на ядра интегрального уравне- 
ния, после чего производится операция обратного БИФ с нормировкой результата. 

Ступень БИФ представлена на рис. 3. 
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Рисунок 3 — Граф-схема одной ступени БПФ 


Комплексные данные для БИФ Ке[тАВ поступают из одного канала. При помощи 
элементов задержки на один операнд 2 раз в четыре такта формируются истинные 
значения КеА, [тА, КеВ, [1тВ. Аналогичным образом получаем коэффициенты для те- 
кущей ступени БИФ Кей’ и [тИ’ из ПЗУ, которое запускается генератором Г`. Резуль- 
таты операции БПИФ формируются каждый четвертый такт. На этом такте в коммутации 
Ксот1 и Ксот2 запишутся реальные части, на первом такте — мнимые части, которые 
специально будут задержаны на один операнд, на втором и третьем тактах — некоррект- 
ные данные. Для объединения двух потоков данных и вырезки из потока некорректных 
данных каждые 2 такта используется блок формирования результата ступени БПИФ 
(БФР БПФ), представляющий собой библиотечный элемент. 

Напомним, что при структурно-процедурной организации вычислений выпол- 
няются все ветви условных операторов. Выбор той или иной ветви осуществляется с 
помощью мультиплексоров (МХ). В частности, мультиплексоры МХ1 и МХ2 выбирают 
ветви с реальными частями только в четвертом такте. В остальные такты коммутиру- 
ются ветви с мнимыми частями. 
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Такая реализация позволяет каскадно подключать ступени БПФ, образуя конвейер. 
Ниже приведен текст программы на языке СОГАМО, описывающий структуру, 
представленную на рис. 3: 


УАБК КетАВ, Кепп\, Ве М, Всоп1, Всопл2, ВеА, 1тА, ВеВ, ПаВ, Ве\, Па\/, ВеАоць ПпАоць ВеВоць ПаВоц : 
геа] сот; 
УАК АВЬиЙ, АВБиР, АВБи 3, АВБо М, УБо А, У/БчЕ2, Пт АБаЕ, НиВЬчЕ : аггау геа] [2:5 геат] сот; 
УАК КК : патьег; 
УАК КК4 : пцеоег сот; 
АВЬий[КК]:=ВевтАВ; 
АВЬиР2[КК]:=АВЬа[КК-1]; 
АВЬиЗ[КК]:=АВЬяЕ?2[КК-1]; 
АВЬ[КК]:=АВЬа{3 [КК-1]; 
ВеА:=АВБо [КК]; тА:=АВЬа [КК]; 
ВеВ:=АВБаР2 [КК]; НиВ:=АВЬа# [КК]; 
МЪоЕ[КК]:=Кепи\/; 
МиР? [КК]:=УБо [КК-1]; 
Ке\!:=\Ъиа2 [КК]; Па\У:=УУБо [КК]; 
ВыЕЕТ(ВеА, ВпА, ВеВ, ПтВ, Ве\\, Пп\/, КеАоцф, ПпАойь, ВеВои, НиВои®); 
ПпАБи{КК]:=ШпАойе; БаВБоКК]:=ПаВойв 
Юг КК:=0 №0 512 4о Берт 

КК4:=КК апа 3; 

ШЕКК4=З Феп Бест 

КВсоп1 :=ВеАой(; Всоп2:=КеВоиб; 
епа; 
е15е Берт 
Всоп1 :=ПпАБа{КК-1]; Всоп12:=ПаВФаКК-1]; 

епа; 
еп4; 
ВешМ:=МакеОпееат(Всот1 ,Ксоп2,512);/БФР БПФ 


Подкадр ВиЕЕТ описывает вычислительную часть структуры БПФ (рис. 3). 

Операция БПФ над трассой выполняется один раз. Поэтому в первом кадре все 
трассы проходят блок вычисления БПФ. После этого результат БИФ каждой трассы 
комплексно умножается на каждое ядро интегрального уравнения, обрабатывается в 
блоке обратного БПФ и нормируется. Количество этих кадров будет определяться числом 
ядер. Следует отметить, что обратное БПИФ можно выполнить на том же аппаратном 
ресурсе, который вычисляет БИФ, просто подав на схему другие коэффициенты. 

Реализация модулей «Вычитание прямого сигнала» и «Обращение интеграль- 
ного уравнения» на РВС с использованием языка СОГАМО показала достаточно вы- 
сокую производительность РВС, которая составила порядка 50 Гфлопс. Выигрыш по 
времени решения обеих подзадач по сравнению с обычным персональным компью- 
тером составил около 100 раз. 
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Розв?язання задач! длагностики дорожнйх покриттв на реконфигурованй 

обчислювальнй систем! 13 застосуванням мови СОГАМО 

У дашй статти розглядаеться розв’язання задач1 дагностики дорожн!х покритив на реконфигурован!й 
обчислювальний систем! 13 застосуванням мови СОГАМО. Детально розглянут! особливост! реал1заци 
обчислювань по трудомпстких частин задач, представлен! граф-схеми, наведен! приклади тексту програми 
мовою СОГАМО. 
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